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ABSTRACT 

We derive a new formalism for convective motions involving two radial flows. This 
formalism provides a framework for convective models that guarantees consistency 
for the chemistry and the energy budget in the flows, allows time-dependence and 
accounts for the interaction of the convective motions with the global contraction or 
expansion of the star. In the one-stream limit the formalism reproduces several existing 
convective models and allows them to treat the chemistry in the flows. We suggest a 
version of the formalism that can be implemented easily in a stellar evolution code. 

We then apply the formalism to convective Urea cores in Chandrasekhar mass 
white dwarfs and compare it to previous studies. We demonstrate that, in degenerate 
matter, nuclear reactions that change the number of electrons strongly influence the 
convective velocities and we show that the net energy budget is sensitive to the mixing. 
We illustrate our model by computing stationary convective cores with Urea nuclei. 
Even a very small mass fraction of Urea nuclei (as little as 10~ 8 ) strongly influences 
the convective velocities. We conclude that the proper modelling of the Urea process 
is essential for determining the ignition conditions for the thermonuclear runaway in 
Chandrasekhar-mass white dwarfs. 
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1 INTRODUCTION 

Gamow & Schonberg (1941) were the first to recognise the 
Urea process (electron captures and emissions on pairs of 
nuclei that can be converted into each other by electron 
captures/beta decays) as a potentially strong source of neu- 
trino cooling in degenerate stars. This process is already 
responsible for significant cooling during the late radiative 
phase of accreting C/O white dwarfs. For each Urea pair, 
the cooling occurs at a mass shell, a so-called Urea shell, de- 
termined by the characteristic density for the pair at which 
the electron captures/beta decays take place. When carbon 
burning starts, a convective core grows and soon engulfs the 
Urea shells. The convective motions across the Urea shells 
back and forth directly affect the net energy release as well 
as the net amount of electrons captured. The resulting phe- 
nomenon is called the convective Urea process and is a key 
ingredient in linking the late evolution of the progenitor of 
a Type la supernova (SN la) with the subsequent explosion 
(Paczyhski 1972). 

Over the last 30 years there have been numerous studies 
of the convective Urea process with often mutually exclusive 
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conclusions. Bruenn (1973) realised that nuclear heating in 
Urea matter outside chemical equilibrium dominates over 
the neutrino losses. On the other hand, Couch & Arnett 
(1975) stressed the cooling effect of the work done by con- 
vection. In the most detailed work to date, Iben (1978a,b, 
1982) computed the evolution of an accreting white dwarf 
including the detailed chemistry of many Urea pairs. He re- 
alised that the turn-over time scales would be of the same 
order as the chemical time scales for the Urea reactions. 
This, he concluded, implied that the mixing processes caused 
by the growth of the convective core would affect the heat- 
ing/cooling by the Urea process. However numerical prob- 
lems caused by his treatment of the mixing prevented him 
from following the computations up to the thermal runaway. 

Barkat & Wheeler (1990) revisited the problem of the 
convective Urea process and provided a clear summary of 
the convective Urea mechanism though later Mochkovitch 
(1996) and Stein, Barkat & Wheeler (1999) pointed out 
a mistake in their treatment and argued that more atten- 
tion needs to be paid to the kinetic energy of convection. 
Bisnovatyi-Kogan (2001) showed that the feedback of the 
Urea process on convection itself should be taken into ac- 
count. As this summary shows, a consistent picture of con- 
vection that properly treats the chemistry is still missing 
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and is badly needed in order to address the problem of the 
convective Urea process. 

Such a theory for chemistry coupled with convection 
was attempted by Eggleton (1983) with the help of a "sim- 
ple rule-of-thumb procedure". Later, Grossman et al. (1993) 
produced a physically consistent model of convection which 
includes chemistry based on a statistical approach for the 
convective blobs. However, their model did not necessarily 
conserve energy and has only been checked without chem- 
istry. 

Here we devise a model for convection that ensures 
energy conservation. We start with the conservation equa- 
tions of radiation hydrodynamics in a spherical configu- 
ration. Using a simple geometry, which mimics convective 
rising plumes, we derive a physically self-consistent model 
of convection which includes time-dependent chemistry. We 
then compare this formalism to previously derived models of 
convection. In the process we obtain a model for the convec- 
tive Urea process that addresses all the problems mentioned 
above, nuclear heating, mixing, convective work, kinetic en- 
ergy and the feedback of the Urea process on convection. To 
illustrate the formalism, we apply it to stationary convective 
Urea regions and show how it affects the energy budget, the 
convective properties and the chemical stratification. In a 
follow-up paper, we plan to apply the formalism in a realis- 
tic, time-dependent stellar model, coupled with a complete 
nuclear reaction network, to determine the ignition condi- 
tions for the thermonuclear runaway in a SN la. 

In section 2 we derive the basic equations for the two- 
stream formalism and suggest a simple model for the ex- 
change of matter, momentum and energy between rising 
and descending flows. In section 3 we compare our model 
to existing models of convection. In section 4 we describe 
the convective Urea process in view of our model. In section 
5 we compute stationary Urea convective cores. We discuss 
and summarise our results in sections 6 and 7. 



2 A TWO-STREAM FORMALISM FOR 
CONVECTION 

Cannon (1993) designed a two-stream algorithm to post- 
process the evolution of chemical species in convective re- 
gions. This model was well suited to the study of convec- 
tive regions in which the chemical time scales were shorter 
than the convective turn-over time scales. But, as a post- 
processing algorithm, it did not tackle the feedback effects 
the chemistry could have on the convection. Here we ex- 
tend his ideas to all state variables in the two streams and 
explore the interactions between mass, energy, momentum 
and chemical transfers between the streams. 

Let us consider a sphere of gas with purely radial ve- 
locities and stratified properties. We assume that, on a shell 
of radius r, there are two different velocities associated with 
the two streams. We further assume that, on this shell, all 
gas parcels moving with a given velocity have the same state, 
they have homogeneous temperature, pressure and chemical 
composition. Figure (1) schematically shows the geometrical 
configuration we have in mind. 

One of the two velocities has to be greater than the 
other, and we refer to the gas moving with this velocity as 
the upward moving gas, even though both velocities could 




Figure 1. Geometrical configuration of the two streams. 

be negative in principle. We write the two velocities as v + 
v u and v — «d, where v u and Vd are both positive and are 
associated with the upward and the downward moving fluid, 
respectively. The ambient velocity v will be defined further 
below. 

We now affix suffixes 'u' to upward moving fluid proper- 
ties and 'd' to downward moving fluid properties. We assume 
that the sound-crossing time of a horizontal section of one 
stream is short compared to any other time scale and hence 
make the approximation that the pressure p — p u — pd is 
the same in the upward and downward moving fluids. The 
state of the gas in both fluids is now completely determined 
by the temperatures T u and Td and abundance vectors N u 
and N d (number of particles per unit mass), provided we 
know the equation of state of the gas. The latter equation 
provides the mass densities p u (p, T u , N u ) and pd(p, Td, Nd), 
and the specific energies e u and e d , along with the radiative 
volumic energies E u and E d - 

We define S u and Sd as the areas of the surfaces occu- 
pied by the two fluids at the shell of radius r. Therefore 

S = S u + S d = 4vrr 2 (1) 

and we define the velocity v by setting the mass flow 

rh = S u v u pu = SdVdPd- (2) 

This defines v as the radial velocity of the centre of a mass 
shell. The net mass flow through the shell moving at velocity 
v is zero. 

We now write the equations for the variation of mass, 
momentum and energy for each of the two fluids and define 
the exchange terms. We then compute the equations for the 
average fluid and the specific equations of evolution for both 
fluids, before computing the differential evolution between 
the two fluids. Finally, we propose a very simple model for 
the exchange terms between the streams. 

2.1 Conservation equations 

We consider the mass, momentum and energy on a shell at 
radius r for each of the two fluids. This allows us to define 
the exchange terms in a conservative way and makes it easier 
to derive the mean equations. Moreover, with this approach 
the specific exchange terms can be defined more rigorously. 
The viscosity of the fluid is neglected as well as the 
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molecular diffusion. All horizontal effects are implicitly in- 
cluded in the exchange terms. 



2.1.1 Mass 

To simplify the derivation, we neglect the mass changes due 
to the nuclear transformations and assume that mass is per- 
fectly conserved. We treat the corresponding nuclear energy 
production only in the energy equation. The rate of change 
of mass is hence equal to the sum of the mass flux in the 
radial direction and sideways, 



^- t (S uP u) = -^{Su P u(v + Vu)]+M, 
and 

d d 

Ol(SdPd) = --^[S d pd(v - v d )] - M, 



(3) 



(4) 



where M represents the mass per unit radius and unit time 
exchanged in the shell from the downward to the upward 
moving fluid. 



2.1.2 Momentum 

The rate of change of momentum is the sum of the momen- 
tum flux (ram and thermal pressure) in the radial direction 
and sideways, added to the gravitational forces where the 
gravitational potential is assumed to be spherical, 

^[S u pu(v+v u )] = -Su^[pu(v+v u ) 2 +p]-S u pu^r+Q,(5) 
and 



dt 



[S d Pd{v-v d )] = ~S d ^[p d (v-v d f +p]-S d p d ^^-Q.(6) 

Here G is Newton's constant of gravitation, m is the total 
mass contained inside the sphere of radius r, 



= I (Supu 

Jo 



+ S d p d ) dr' 



(7) 



and Q is the momentum per unit mass and unit time ex- 
changed in this shell from the upward to the downward mov- 
ing fluid. 



2.1.3 Energy 

The rate of change of energy is the advective plus diffusive 
energy flux, in the radial direction and sideways, added to 
the work of pressure forces and heat generation, 

^(SuPue u ) = -^[S u p u e u (v + v u )] - p-^-[Su(v + v u )] 



-p-Ki(Su) + S u pu I e u - 

dV ' V S u p u dr 



+ E 



(8) 



and 

^(S d p d e d ) = -JL[S dPd e d (v-Vd)]-p-l^[Sd(v-Vd)] 



d , . / dL d 

-p-^TWd) + S d p d e d ~ — — 

dV \ SdPddr 



Here, L is the luminosity carried by the conductive and ra- 
diative processes (in the radial direction), 



and 
L d 



S u c dE u 
3k u dr 

PuCPuXu-T^Ou, 



_S d cdE d 
3n d dr 

0T dQ 

PdCpdXd-7^7<3d, 



(10) 
(11) 

(12) 
(13) 



where E is the radiative energy per unit of volume given 
by the equations of state, c is the speed of light, k is the 
flux weighted total opacity for the conductive and radiative 
processes (in the diffusion approximation), \ is the thermal 
diffusion coefficient and cp is the heat capacity at constant 
pressure. 

The net energy production per unit time and per unit mass 
is e. This can include the nuclear energy production if the 
rest mass energy is not yet included in the specific energy, 
and it includes the neutrino losses. The energy per unit time 
and per unit radius exchanged in the shell from the upward 
to the downward moving fluid is E. 



2.1.4 Chemistry 

The rate of change of a species is its flux in the radial direc- 
tion and sideways added to its chemical rate of change, 

^(5* uPu N u ) = ~-?-{S u p u N u (v + v u )] + S u p u TL u + -N (14) 

and 

d ,„ „ , d 



dt 



(SdPdNd) = - — [S d p d Nd(« - v d )] + SdPdRd - N, (15) 



where Nj represents the mass of species j per unit radius 
and per unit time exchanged in the shell from the downward 
moving fluid to the upward moving fluid, and Rj is the rate 
of change of species j per unit time and unit mass due to 
nuclear reactions. 



2.2 Equations for the mean fluid 

When we sum up each pair of the above equations we obtain 
the equations for the variation of the mean mass, energy and 
momentum. We first derive these equations in their volumic 
form and then gather them all in a more familiar specific 
form. 



2.2.1 Mass 

The mean density p is defined by 

pS = p u Su + p d S d . (16) 
The equation for mass conservation then reads 

(17) 



§ t (Sp) = -l(Spv), 



E. 



(9) 



which is the unchanged continuity equation for the mean 
fluid. 
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2.2.2 Momentum 

The mean equation for momentum is more complex and 
involves several additional terms, 



-(Spv) = S—(pv 2 + p) - Sp-^ + SpX c . 



dt 
where 



dr 



S pX conv = X\ + X2 + Xz 

x l 



Su^{puvl) - S d -^(pavj), 



d d 
X 2 = -2v[S u -^{p u v u ) - S d -^{p d v d )] and 



X a = (p d - Pu)v 2 



(18) 

(19) 
(20) 

(21) 
(22) 



S u 9S d Sd dS u 

where we have used equation (1). Usually, the evolution of 
the star is quasi-static, which means that v is negligible. 
Then only X\ should be retained because X2 and X3 are of 
order 1 and 2 in v. 

Xi accounts for a convective pressure support in the 
flow. It slightly changes the hydrostatic equilibrium but is 
usually negligible in the subsonic regime. 

2.2.3 Energy 

We define the mean specific energy as 

Spe = p u S u e u + PdSded (23) 

and the mean energy generation per unit mass e in the same 
way as 

pe = p u S u e u + pdSd^d- (24) 

The radiative and conductive luminosity is given by 

L = L n + L d . (25) 

Then the equation for the mean energy is rather simple and 
includes only two additional terms, 

§- t (Spe) = -l(Spe V )- P l { S V ) 



+Sp 



where 

Lconv 



(.- 



dL + L co 
Sp~dr~ 



= m(h u — h d ) = -Spu(h u — /i d ), 



(26) 

(27) 
(28) 



VVconv = -U( )— . 

2 p u pd or 

Here h — e + p/p is the enthalpy and u = 2rh/(Sp) is the 
mean convective velocity. 

Note that the convective luminosity L conv naturally appears 
as an enthalpy flux. 

The term W™ can be interpreted as the work done by the 
buoyancy forces. Because the pressure is decreasing upward, 
and usually p u < pd, it always provides a sink term for the 
energy. 

Note that this approach cannot be used to write down 
a general equation for the entropy because the thermody- 
namic relations are only valid for the up and down streams 
separately but not for the mean fluid. We will derive an en- 
tropy equation later in the one-flow limit (see section 4.3). 



2.2.4 Chemistry 

The equation for the mean abundances is 

Wt {SpN) = -|^ Nv ) + SpK - ^^ Fconv ' (29) 
where 

SpN = S uPu N u + S d p d N d (30) 
and 

F conv = m(N u -N d ), (31) 
which we interpret as a diffusion flow in section 2.5. 

2.2.5 Mean specific equations 

We now summarise all of the mean equations derived in the 
previous section in their specific form (i.e. per unit mass), 
where we use the notation 2| = §| + «ff and div(x) = 

1 dSx . 
S dr ■ 



1 Dp 
p Dt 

Dv 
Dt ~ 



— — div(v) 



2^ 



1 dp Gm 
p dr r 2 



+ X C , 



De D 1 _ dL 
Di +P W~p~> ~ £ ~ dm~ 

and 

DN_ d 

-Dt- K ~dm~ hconv 

with 



8L C 



dm 



+ W a 



dr 



X C onv 

Lconv = m(h u -hd), 



1 

S~p 



Su-^-(puvl) + Sd-^-(pdvj) 



dr 



W COI 

and 
F 

- 1 - con - 



1 

= 2 U 



dp 

Pu Pd J dr' 



m(N u -N d ). 



(32) 
(33) 
(34) 

(35) 

(36) 
(37) 
(38) 

(39) 



These equations are the usual equations of radiative 
stellar evolution with additional terms due to the differen- 
tial motions in the two fluids. The latter terms depend on 
the convective velocities v u and v d , as well as on differences 
between quantities in the two flows. We derive equations for 
the mean convective velocities and those differences in the 
one flow limit in section 2.4. 



2.3 Specific equations 

We first rewrite the conservation equations in their specific 
form. This will allow the derivation of the equations in the 
one-flow limit when we take the difference between the spe- 
cific equations for the up and down motions. 



2. 3. 1 Momentum 
8_ 
dt 



d 

(V + V u ) + (v + ll„) — (l) + V u ) - (V + V u ) 



1 dS u 
Su dr 
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_J_^_Gm + l ^_ Kt 24A Mass 

p u Or r 2 S u pu 

Taking the difference between equations (3) and (4) and 
and using equation (32) we find 

|(« - t*) + (v - v^v - „ d ) - („ - , d ) 2 i-fl = D = 16 _ 2M 



1 dp Gm 1 . • ■ . ,. . . 

~p~~dr r 2 S p ~ ~ This equation describes the evolution of the asymmetry An 



Dt \ u ) Spdr y ' ' Sp 

This equation describes the evoh 
of the drift motions relative to the mean velocity v. 



2.3.2 Energy 
<9e u de u 



d 1 . ,91 

^7 h (V + V u )- 

at p u Or p u 



Supu dr S u pu 
and 

de d de d 
_ + (u _ Ud) _ +p 



(42) 



d 1 , ,91 

^7 h ( v ~ Vd >7i 

at pd or p d 



1 dLd l_ 

SdPd dr SdPd 



(E-Mhd) 



(43) 



2.3.3 Chemistry 

^ + (, + , u )^R u + ^(N-MN u ) 



and 

9Nd 

dt 



2.4 



+ (^« d )^ = R d - tt^(N - MN d ). 



dr 



SdPd 



(44) 



(4£ 



The difference equations in the one stream 
limit 



In the following we make the approximation that the relative 
difference of a quantity between the two streams is small. For 
each quantity x we write x u = x + Ax and Xd — x — Ax 
with Aa; <§C x. We note that the previously defined mean 
convective velocity u is in fact the arithmetic mean of v u 
and Vd ■ Hence, we write v u = u + Au and Vd = u — Au. 
In contrast, we keep the former definition for S as the total 
surface of the shell of radius r so that = \S + AS and 
Sd = t\S — AS. We then compute the difference of the spe- 
cific equations (40)-(45) between the two fluids and neglect 
second order terms in the A quantities. 

In this approximation, the averages defined in section 
2.2 above are all arithmetic means. We also note that the 
standard thermodynamic relations hold for the A quantities, 
which greatly helps when deriving the entropy equation (see 
section 4.3). 

Equation (2) yields 



AS Ap 
S p 



Au 



= 



(46) 



and this allows us to eliminate AS in what follows. It can be 
thought of as the equation governing the convective motions 
with u as the convective velocity. 



2.4.2 Momentum 

The difference between equations (40) and (41) yields 



Dm / d(v + 2Au) v + Au 
Di +U { dr ~ 4 ^— 



Ap dp 

T 2 



d_ 

' dr 



Ap 
P 



d (Au 



dr 



M(v + 2Au)]. 



(48) 



This equation describes the time evolution of the convective 
velocity u. The main source term is the acceleration due to 
the buoyancy force while the main sink term comes from the 
momentum exchange term between the two flows. 
Apart from the Au terms, — (v 2 +u 2 )-^{^) is the only non- 
local term. That is it is the only term that involves spatial 
derivatives of convective properties. 



2.4-3 Energy 

The difference between equations (42) and (43) yields 

DAe D / Ap\ ( de 8 l\ 
^~ P Vt[-p 2 -) +U \dr +P d-r-p) = 



Au dL dAL 2 , ■ 

Ae + 2 — + 2—— + —{E - Mh). 

u dm dm bp 



(49) 



This equation implicitly determines the time evolution of 
the temperature difference between the two flows as we show 
in section (3.1). Note that the At term may help to either 
increase or decrease such a temperature difference, depend- 
ing on the temperature dependence of the energy genera- 
tion rate. The other non-local term in the energy equation, 
2^T' accounts for differenital thermal diffusion along the 
two columns of fluid. 

2.4.4 Chemistry 

The difference between equations (14) and (15) yields 



DAN d~N ._ 2 ,-r. 

— + ^=AR+^(N-MN). 



(50) 



This equation describes the evolution of the chemical compo- 
sition difference between the two flows. In some situations, 
the temperature dependence of the nuclear reaction rates 
may affect this difference through the AR term. 
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2.5 A model for the exchange terms 

All horizontal motions and transport phenomena are mod- 
elled through the exchange terms. These need to be specified 
to close the two systems of equations (32)-(35) and (47)- 
(50). With these, the equations provide a complete time- 
dependent model for convection which guarantees the con- 
servation of mass, momentum, energy and chemical trans- 
formations. Most of the existing convective models provide 
approximations for these exchange terms. Here we present 
a simple, somewhat ad hoc, but physically plausible choice 
for these exchange terms, 

Q = -{v u + v d )rh/\ + Mv, 



E = (/la - h u )m/X + , 
and 



Cp(T d 



7YI 

T u )-+Mh, 



N = (N d - N u )m/A + MN. 



(51) 
(52) 



(53) 



The first terms in these expressions account for momentum, 
energy and chemical exchanges without any net transfer of 
mass. They are designed to mix the two fluids on a length 
scale A. The second term in the energy exchange equation 
(52) accounts for horizontal heat diffusion across the edges 
of the streams. The parameter /3 is a form factor which fixes 
the ratio of the perimeter of the streams to their separation 
times S/X 2 . We use [3 — 9/2 to recover the exact formula- 
tion of classical mixing-length theory (MLT; Bohm-Vitense 
(1958) as presented in Kippenhahn & Weigert (1990)). The 
last terms in expressions (51)-(53) are the fluxes due to a 
net transfer of mass M from one stream to the other. Such 
transfer must exist at the outer boundaries of a convective 
zone where the fluid effectively makes a U-turn. 

There is no obvious, simple physical prescription for M 
so we make the assumption that u — v u = t> d and convective 
motions are symmetric with respect to the mean velocity v. 
This approximation is in fact implicit in almost all convec- 
tive models to date. The term M is then given by equations 
(2)-(4) as 



M 



dm 
dr ' 



(54) 



We now have a complete time-dependent description of the 
convective properties of the flow and can investigate the 
characteristics of this model. 



3 COMPARISON WITH EXISTING MODELS 
FOR CONVECTION 

In this section we consider the approximations made in a 
number of theories of convection. Using the same assump- 
tions, we derive the equations for the convective motions in 
our framework and emphasise the characteristics that are 
peculiar to our formalism. 



3.1 Mixing length theory 

Mixing length theory (MLT) assumes a stationary state for 
convection and makes the quasi-static approximation, v ~ 0. 
In this case, all ^ terms can be set equal to zero in the 
difference equations. In addition, we set Aw, Ae, AL and 



AR to zero because the processes that lead to these terms 
are usually neglected in MLTs. We also neglect the non-local 
term in the momentum difference equation. The difference 
equations then simplify to 

„,2 



Ap dp 
p 2 dr 
de + d_l 
dr ^ dr p 
dN 
dr 



Aft vepAT 
" 2 - " and 



= -2 



AN 



(55) 
(56) 
(57) 



Note that equation (55) is not exactly the same as in 
classical MLT : there is usually a factor of 1/2 multiplying 
the left hand side of this equation to account for the fact 
that half of the work done by the buoyancy force is used 
to push aside the surrounding medium when a convective 
element rises. 

From these equations we extract the convective velocity, 
the temperature and chemical differences between the two 
streams (see appendix Al.l for a detailed derivation). These 
can now be used to express the convective terms in the mean 
equations, 

1 d 



X, 



p dr 



(pu ), 



1 



m = -Spu, 



m{c p AT + fi'.AN), 

3 



conv 

F conv = 2rhAN 



^— and 
A 



(58) 

(59) 
(60) 

(61) 
(62) 



where fi' — T(-^) T , P + (J- and /x is the chemical potential 
(see appendix Al.l). 

With these equations, we have derived a MLT that 
is consistent with chemistry. While the excess temperature 
equation is unchanged, the convective velocity now depends 
on the chemical stratification (through the cubic A5) and 
so does the convective luminosity. There is also an addi- 
tional work term due to the fact that we assume a reversible 
process for the momentum exchange (see section 3.3). We 
also obtain an explicit change in the condition of hydro- 
static equilibrium owing to convection. Finally, convection 
naturally appears as a diffusion process for the chemistry. 

3.2 Unno (1967) 

If we include the D/Dt terms in equations (55) and (56) 
and further neglect the chemistry dependence of the con- 
vective luminosity and velocity we immediately recover the 
same time-dependent version of MLT as Unno (1967) for the 
excess temperature and the convective velocity. 



3.3 Kuhfuss (1986) 

Kuhfuss (1986) only computes the evolution equation for the 
convective velocity. He uses a diffusion model to compute the 
correlations between velocity perturbations and any other 
perturbation. This model is recovered in our formalism if 
we set ^ = and AR = Ae = AL — in the difference 
equations for energy and chemistry. 

Furthermore, each term of his equation (25c) for the 
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convective kinetic energy corresponds to one term in our 
equation (48), except that his non-local term ^^divj t is 
different from ours (see paragraph 3.5), and we do not ac- 
count for his viscous terms (we assumed an inviscid fluid). 

Finally, if we compare our equation for the average in- 
ternal energy of the gas, we note that our model misses the 
heat production owing to dissipation of convective motions. 
This is due to the fact that we assume a reversible exchange 
of momentum so that there is no associated heat production. 
In contrast, Kuhfuss (1986) assumes viscous dissipation for 
the convective motions. 



3.4 Eggleton (1983) 

Eggleton (1983) uses his rule of thumb to average the hydro- 
dynamical equations in order to get evolution equations for 
the mean fluid and the perturbed quantities. He then obtains 
a full set of equations that can be identified with our mean 
and difference equations. Our formalism agrees fairly well 
with this rule of thumb and hence provides a more physical 
basis for it. 

We can recover almost the same equations if we set 
Am = 0. The local models differ only in the mean energy 
equation because Eggleton neglects the thermal part of the 
chemical potentials. However, those terms are important 
when beta decays or electron captures occur in degenerate 
matter. 

As for the non-local terms, in the energy difference 
equation takes the same form of the term in equation (41) 
of Eggleton (1983) that accounts for the thermal diffusion 
through the front and back of the eddies. However, our non- 
local term in the momentum difference equation does not 
agree with equation (42) of Eggleton (1983) for the velocity 
perturbation. 

An interesting point is that Eggleton (1983) found a 
term similar to our AR term. We have shown that an addi- 
tional term Ae enters into the energy equation. These terms 
account for differential reactivity in the two streams. How- 
ever, our stationary computations show that these effects are 
negligible as far as the convective Urea process is concerned 
(see section 6). 



3.5 Grossman et al. (1993) 

Grossman et al. (1993) use the Boltzmann equation cou- 
pled with dynamical equations to compute a hierarchy of 
moments for the hydrodynamical equations. In their frame- 
work, our formalism can be recovered if we specify the dis- 
tribution function of the blobs as 

f A (t,z,v,T) = ^[S u p u 5(v+v u ,T u )+S d p d 5(v-v d ,T d )]. (63) 

This allows a direct comparison between the two formalisms. 
With this definition, our variables v, p and T correspond to 
their variables v, p, and T. 

The zeroth-order equations of their hierarchy can be 
directly compared with our mean equations. The velocity 
equation is found to differ by terms negligible in the subsonic 
regime. Their temperature equation without chemistry has 
an additional source term due to viscous dissipation just as 
Kuhfuss (1986). Their temperature equation with chemistry 
misses correction terms in the luminosity, the work term 



and the source term. This is due to their use of a dynami- 
cal equation for the entropy which does not account for the 
changes owing to chemical evolution. This suggests that en- 
ergy conservation may not hold in their case when chemistry 
is included. 

Their higher-order equations can also be compared with 
our difference equations, although this is less straightfor- 
ward. We did such a comparison but only for the velocity 
difference equation. We then obtain the same non-local term 
as Kuhfuss (1986): their 4J|(pw 3 ) corresponds to ^^divjt 
of Kuhfuss (1986). If we substitute our distribution function 
/a in this term, we obtain 



--(pu An) 

which should be compared to our 

2 dAu 3 d Ap 

2u ~5 hM ^ 

or or p 



(64) 



(65) 



when v = 0. Our term contains Ap and theirs contains the 
spatial derivative of the convective velocity. On the other 
hand, our treatments agree as far the chemical dependence 
of the convective velocity is concerned. 



4 THE CONVECTIVE URCA PROCESS 

In this section we first present the basic nuclear reactions 
responsible for the Urea process. Because previous studies 
have mainly concentrated on the entropy equation we derive 
it in our formalism and then examine the influence of the 
Urea nuclei on convection. 



4.1 Urea reactions 

Urea reactions involve pairs of nuclei of the form (^ +1 MjD) 
where A is an odd number. The member of a pair with an 
additional proton is called the mother (M), while the other 
one is referred to as the daughter (D). Electron capture and 
beta decay turn one into the other : 



Electron capture : M + e — > D + v 
Beta decay : D — > M + e~ + v 



(66) 
(67) 



Tsuruta & Cameron (1970) give the reaction rates for 
electron capture and beta decay A + and A - per nucleus. 
They also provide the corresponding neutrino losses L + and 
per nucleus. The typical time-scale for Urea reactions is 
10 5 s. These reaction rates depend mainly on the chemical 
potential p e of the electrons. This is mainly on the mass 
density p of the degenerate matter. Each Urea pair has a 
threshold energy /i t h above which significant electron cap- 
tures can occur. We approximate the Coulomb corrections 
to these threshold energies in the same way as Gutierrez et 
al. (1996). When p e > p t ^ (i.e. p > p t h) electron captures 
quickly turn the Urea matter into daughter nuclei. When 
Pe < pth (i.e. p < p t h) beta decays quickly turn it into 
mother nuclei. When p e ~ p th (p ~ p th ) both reactions are 
significant and the Urea matter quickly evolves into a mix- 
ture of mother and daughter nuclei. A shell on which p = p t h 
is called an Urea shell. 

When the Urea matter is in chemical equilibrium both 
reactions take place at the same rate, which is highest near 
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the location of the Urea shell. Because both reactions emit 
neutrinos, this leads to neutrino cooling which is strongest 
at the Urea shell. When the Urea matter is far from equilib- 
rium, nuclear heating takes place at the same time and usu- 
ally dominates over the neutrino cooling (see Bruenn 1973, 
and section 4.3.1 below). 

The pair 23 Na/ 23 Ne is a typical example and perhaps 
the most active Urea pair in massive white dwarfs. Its 
threshold energy is /i t h = 4.38 MeV, which corresponds to 
pth = 1.7 x 10 9 g.cm~ 3 in purely degenerate matter. This 
density is slightly below the density for carbon ignition in 
a C+O white dwarf accreting at a rate of 10~ 7 MQyr -1 . 
As a consequence, the growing convective core soon engulfs 
the corresponding Urea shell after carbon has ignited in the 
centre. 

Here we assume that the stellar matter contains a fixed 
number Ntj of Urea nuclei (per unit mass) of a given pair. We 
define Nm and TVd as the corresponding number of mother 
and daughter nuclei so that iVu = Nm + No- Then the 
rate of change of mother and daughter nuclei per unit mass 
becomes 



-Rd = -X + N M + X~N U . 



(68) 



If we write N M = N v A /(A+ + A ) = N v - 7V D as the num- 
ber of mother nuclei per unit mass at chemical equilibrium 
we can rewrite the reaction rate as 



Rb 



(69) 



Useful explicit expressions for the various thermody- 
namical properties of the Urea nuclei are 



fiM = (Z + l)/i e + fcT(a + lniV M ), 

/id = Zfi e + kT(a + In 7V D ), 

fl' M ^ (Z + l)/ie + i*T, 

Md ~ Zfi e + X -kT, 

Mm = -(Z + l)Ar M /iV e and 

Md = -ZN D /N e , 



(70) 
(71) 

(72) 

(73) 

(74) 
(75) 



where N e is the total number of electrons per unit mass, k 
is the Boltzmann constant and a is a combination of tem- 
perature and density logarithms (see Pols et al. 1995). The 
approximations for /i M and /t D are obtained by neglect of the 
electron and radiation contributions to the entropy. In the 
strongly degenerate case the kT terms are usually negligible. 

4.2 Convective velocity 

In appendix A (equation A8) we derive the approximation 
which relates the convective velocity to the temperature and 
chemical gradients 



u = ui y/S{V -V a )- (76) 
where the term depending on the Urea process is 



m".V n = - — [(Z + l)iV M V M + ZJVdVd]. 

iV e 



(77) 



Since JVu is assumed to be uniform, ATmVm + -/VdVd = 0. 
Hence relation (77) becomes 



,," V7 - Nm T7 
fjL .Vn = — Vm- 



(78) 



Electron captures are much stronger in the centre of the star 
than in the outer regions so the mother fraction generally 
increases outwards, /j".Vn > and the effect of the pres- 
ence of Urea pairs is to reduce the convective velocity. In 
sections 4.4 and 5 we show that this effect can actually be 
quite strong and may even inhibit convection. Note that, in 
the case where Ny is not uniform, the sign of the /j".Vn 
term may change. 

4.3 Entropy equation 

One of the biggest uncertainties in previous models for the 
convective Urea process was the form of the equation for 
the evolution of the entropy s in the presence of convection. 
Here we derive the entropy equation in the one-flow limit. 
The small A approximation allows us to use thermodynamic 
relations for the mean fluid and we can transform the left- 
hand side of equation (34) so that 

dL dL c 
dm 



^Bs DN 

T h u, 

Bt ^ Bt 



dm 



(79) 



We then use equation (35) and the relation Ah — /a. AN = 
TAs (because Ap = 0) to write the equation for the mean 
entropy as 



T 



Ds _ , dL 
Bt ~ £ dm 
where we define 



dL' 



dm 



+ W' 



e = e - £t.R, 

L' = 2mTAs 



and 
W' 

' ' cor 



2mAN 



dfi 

dm 



(80) 

(81) 
(82) 

(83) 



The entropy equation then takes a form that is similar to 
the energy equation but with different definitions for the net 
heating, the convective luminosity and the work. We next 
consider how these terms are affected by the chemical state 
of the Urea matter. 

4-3.1 Net heat generation 

The net heating due to Urea reactions is 



ev = Mth-RM — L Nm — L No - 

In the very degenerate case, /iM 
we obtain 

e' u = N v C+{N M -N M _)H, 



^mRm — hdRd- (84) 
^ (Z+ l)/i e and /id ^ Zfi e , 

(85) 



where 



C = - 



L+\- +L-\+ 



H = -L 



A+ + A- 

' + L~ + (/i e - /ith)(A" t 



+ A-). 



(86) 
(87) 



In figure 2 we plot \C\ and \H\ for the Urea pair 23 Na/ 23 Ne. 
The net Urea heating t' v is the sum of two terms. At chem- 
ical equilibrium, only the first term remains. This is always 
negative and so causes cooling. The sign of the second term 
depends on the signs of A^m — N M and H. The term Nm — N M 
is likely to be positive below the Urea shell and negative 
above it. The value of H has the same sign except close 
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Figure 2. Functions \H\ (solid and dotted lines) and — C (dashed 
line) for the 23 Na/ 23 Ne Urea pair at a temperature of 3 X 10 8 K. 
The solid portions of \H\ indicate where H and Nyi — have 
the same sign (heating) and dotted otherwise (cooling). 

to the Urea shell. Hence, the second term (Am — N^)H is 
positive, implying heating, except close to the Urea shell. 

The relative magnitude of the second to the first term 
is proportional to the departure from chemical equilibrium 
(Am — A m )/Atj. One can see from figure (2) that the heating 
can be balanced by neutrino losses only if the system is 
close to chemical equilibrium because, generally, \H\ >• |C|. 
In a convective region mixing puts the Urea abundances 
slightly (or strongly if convection is very efficient) outside 
equilibrium. The amount of convective mixing is therefore 
crucial for computing the net Urea heating. Generally, the 
net effect is cooling close to the Urea shell and heating far 
away from it. 

4-3.2 Convective luminosity 

The contribution to L' conv associated with the Urea pairs is 

Ly = 2m(/j - fi').AN = 2rhkTln(AD/A M )AAM. (88) 

The quantity one needs to compare with e' v and W c ' onv is in 
fact the mass derivative of the luminosity JL^- . This term is 
usually negligible in the very degenerate case. However, at 
convective boundaries, the derivative |^ can be large and 
make this term relatively more important. 

4.3.3 Work term 

Here we compute the contribution of the Urea pairs to the 
work term in the entropy equation, 

W{j = 2riiAN.^ = 2rhANM -tt— ■ (89) 
dm dm 

Because the Urea matter is richer in mother nuclei in the 
outer parts of the star, AAm > 0, and because the density 
is decreasing outwards, fj, e is decreasing outwards. Hence, 
W{j is always negative. 

We note that W[j = is identical to the work term 

Iben (1978b) designed for his computations although he uses 



it in the energy equation and does not specify what convec- 
tive velocity or luminosity he adopts. 

Finally, no heat production is associated with our chem- 
ical exchange model. Indeed, we assume that a reversible 
process is responsible for the mixing of chemical species be- 
tween both streams. An irreversible process for this chemical 
mixing could give rise to a heating term which could balance 
part (or all) of this additional work term. 



4.4 Criterion for convection at the centre of the 
star 

Let us now assume that we know the composition, tempera- 
ture and density at the centre of a star. We can then derive a 
criterion for whether there is stationary convection and de- 
duce an upper bound to the amount of mixing at the centre 
of a convective Urea core. 

We consider a very small sphere of mass m at the centre 
of the star. In a stationary state the total luminosity at 
the edge of this sphere must balance the energy production 
inside the sphere. If we assume that all the energy is carried 
out by convection we can write the convective luminosity 

em = L conv = 2Ahm = 2(c p TAlnT + fi' .AN)m. (90) 

We obtain a similar expression for the net number of parti- 
cles of each kind flowing away from this sphere 

Rm = F C onv = 2ANm. (91) 

Combining these two equations we obtain the density 
difference at the edge of the sphere 

Alnp = JUL Ls e -M'- R + y . (92 ) 

H 2m y cpT 4^ N 3 J y ' 

A stationary convective state exists at the centre if and 
only if Ap < 0. This translates into an upper limit for the 
mother fraction at the centre of an Urea convective core. 
Using the relations fi'j ~ Zjfj, e and p" = —ZjNj/N e , where 
A e is the number of electrons per unit mass, we obtain 

N ™- N ™ <6 £t{^W 1 -B) 1 =^ = ^.(93) 

where e* is the net heating at Urea chemical equilibrium 
(equal to the heating from carbon burning at the centre) and 
Al and Xl are defined by equation (93) . Al is the minimum 
mass fraction of Urea pairs which has a significant effect on 
convection. It generally has a fairly small value (see figure 
3). 

When Atj 3> Al the Urea pairs must be close to chem- 
ical equilibrium. In other words, the convective core can be 
only very slightly mixed. But if they were in chemical equi- 
librium carbon burning alone would produce a large buoy- 
ancy and hence drive strong mixing. Therefore, the Urea 
composition has to adjust itself to balance the heat from 
carbon burning and produce almost zero convective veloci- 
ties. Hence, the inequality (93) is nearly an equality in prac- 
tice and Am — A M ~ Al, a result that is verified in our 
simulations of stationary convective cores (see Section 5.3). 
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Figure 3. Xl (solid line) and p (dotted line) against fi e — fi t ^ 
for the 23 Na/ 23 Ne Urea pair at a temperature of 3 X 10 s K. 



5 STATIONARY CONVECTIVE URCA CORES 

If we set all Jj derivatives equal to zero in equations (3)- 
(15), we obtain a system of coupled ordinary differential 
equations. If we specify the state variables at the centre of 
the star, we can integrate these equations outward, using 
a shooting method, and calculate the hydrostatic profile of 
the star. 

To obtain a guess for the state variables at the centre, 
we run a very simple time-dependent model of an accret- 
ing white dwarf without Urea nuclei. We then compute the 
stationary state of stationary convective cores for different 
Urea compositions and different versions of our convective 
model. 



5.1 Time-dependent model 

We used the Eggleton (1971) stellar evolution code to calcu- 
late the evolution of a white dwarf, composed entirely of 12 C, 
16 and 20 Ne (with mass fractions 12 C = 0.25, ie O = 0.73, 
20 Ne = 0.02) accreting matter at a rate of 10~ 7 MQ/yr. 
The initial mass was taken as 1 Mq. Only the carbon- 
burning reaction 12 C( 12 C,a) 20 Ne immediately followed by 
12 C(a,7) 16 was taken into account. The convective model 
used in the Eggleton code is standard MLT and we use the 
approximate equation of state of Pols et al. (1995) assuming 
complete ionisation. 

We stop the computation during the carbon flash when 
the convective core has reached a mass of 0.4 Mq. We plot 
the convective velocity of this core against mass in figure 
(4). At this point, the threshold density for the 23 Na/ 23 Ne 
Urea pair is in the middle of the convective core. The central 
density is 2.6 x 10 g.cm and the central temperature is 
3.1 x 10 8 K. 

If we were to integrate the equations of hydrostatic equi- 
librium from this central state using MLT, we would obtain 
a fully convective star. The reason is that the core is actually 
being heated and that the term T^f is non-zero. This term 
is negligible in the carbon-burning region but is significant 
in the outer part of the convective region. However, the ef- 
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Figure 4. Convective velocity profile through the core. The solid 
line is for the time-dependent simulation and the dashed line is 
for the hydrostatic simulation. We also plot in the timc- 

dependent simulation (dotted line). 



feet of convection is to homogenise the entropy profile and 
so T§§ is rather uniform in the convective region. We can 
therefore use its value to offset the nuclear heating when we 
compute the hydrostatic profile. This brings the stationary 
convective profile very close to the time-dependent one (see 
figure 4). 



5.2 The shooting method 

Setting all the terms to zero in the equations for the two 
streams, we obtain a system of coupled ordinary differential 
equations, which we integrate numerically from the centre 
of the star. As central boundary conditions for temperature, 
density, C, O and 20 Ne mass fractions we take the results of 
the time-dependent model. For each value of N\j, we make 
an initial guess for the number density Nm- The estimate 
for Nm is then iteratively improved by successive outward 
integrations until the condition ANm = is satisfied at the 
outer edge of the convective region. 



5.3 Results 

We compute stationary convective cores for different mass 
fractions of Urea pairs X v = 23 x N v = , 10~ 12 , 10~ 9 , 
10~ 8 , 10~ 6 , 10~ 3 . A value of X v = 10~ 3 would require a 
very efficient conversion of 20 Ne into Urea pairs and hence 
gives a reasonable upper limit for the possible abundance of 
Urea pairs. 



5.3.1 Without Urea nuclei 

For a pure C+O+ 20 Ne mixture, we compare the velocity 
profiles of the stationary convective core given by MLT and 
the two-stream model (TSM). In the TSM, we first set the 
work term W com equal to 0. Figure (5) shows that the ve- 
locities differ by a factor of \[2. This is the only difference 
between the two models and it can be traced back to the fac- 
tor 2 in equation (55). Indeed, MLT has a factor 4 instead 
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Figure 5. Convective velocity profile of different models of sta- 
tionary cores without Urea nuclei. Solid and dashed lines are for 
MLT and TSM models with viscous dissipation of momentum. 
Dotted and dash-dotted lines for MLT and TSM models with the 
work term W C onv 



Figure 6. Fraction of mother nuclei Xm/^v against mass in 
the upward and downward moving streams for X\j < 10~ 9 . The 
mean is also indicated. The model is the two-stream model (TSM) 
including the work term. 



because it assumes that half of the work done by the buoy- 
ancy forces is used to push aside the surrounding medium 
when a convective element rises. 

When we put the work term back in MLT or TSM the 
resulting convective core shrinks. Less convection is needed 
to carry out the C burning energy. This is not surprising 
because adding the work term is equivalent to suppressing 
the viscous heat produced by the dissipation of the drift 
motions. 

5.3.2 Very low Urea abundance 

For X\j < 10 -9 , the Urea nuclei do not have an impact on 
the convective velocity. They are mixed passively through 
the convective region. Figure 6 shows the relative abundance 
of mother nuclei Nm/N\j in both streams. Mother nuclei 
come from above the Urea shell. They capture electrons as 
they descend below the Urea shell and are converted into 
daughter nuclei. As these rise back above the Urea shell, 
they emit electrons and the number of daughter nuclei rises 
again. Finally, they cycle back down through the Urea shell. 

The maximum relative difference of composition be- 
tween both fluids is 17%. This suggests that the one-stream 
approximation may well be adequate. 

5.3.3 High Urea abundance 

For Xu > 1CF 8 ~ Xl our shooting method is not able to 
find a stationary state. When the mother mass fraction at 
the centre is too high the relative composition difference 
diverges as the convective velocity tends to zero. When it 
is too low the convective velocity reaches a minimum and 
then rises again up to the edge of the white dwarf. Figure 7 
shows the convective velocity profiles in the latter case and 
demonstrates how drastically the Urea nuclei can affect the 
velocity profile. 

We illustrate this effect by plotting the density differ- 
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Figure 7. Convective velocity profile of different models of sta- 
tionary cores with Urea nuclei. The different values for X\j (the 
total mass fraction of Urea nuclei) are indicated next to each 
curve. 



ence A In p in figure 8 for a case where Xm at the centre is 
too low. Indeed, the density difference controls the buoyancy 
force and hence the convective velocity. This figure shows the 
competing dependence of A In p on the temperature and on 
the chemistry (the Urea nuclei). As was previously noted 
in section 4.2, the chemical part of the density difference is 
positive and stabilises convective motions. In the outer parts 
of the convective region the temperature and chemical de- 
pendence cancel each other to yield a very small convective 
velocity. This suggests that the criterion for semi-convection 
might be fulfilled in the region just above the point where u 
vanishes. 

In degenerate matter the mass density is tied to the 
electron abundance. Hence, all reactions that change the 
number of electrons have an effect on the buoyancy. The 
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Figure 8. Density difference Alnp through the convective core 
for a low central mother fraction case. Alnp = — 5AlnT + 
/i".AlnN. It is the sum of a temperature-dependent term 
and a chemistry-dependent term. We plot — 5AlnT (solid line), 
jU".AlnN (dashed line) and Alnp (dotted line). 
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Figure 9. Quantities e' v (solid), (dotted) and (dashed) 

are plotted for a low central mother fraction case of model X\j = 
10~ 3 (TSM with work term). The total net energy generation 
rate is greater than 5 X 10 4 crg.g _1 .s _1 . It is dominated by the 
energy generation due to carbon burning. 



terms e\j, and W{j seem to be of secondary importance 
compared to the change in the convective velocity. We plot 
them in figure 9. They are relatively unimportant compared 
to the energy generation rate due to carbon burning. 



6 DISCUSSION 

In principle the formalism we have derived is sufficiently gen- 
eral that it can be used with any model for the exchange (or 
diffusion) between the two streams. It automatically guar- 
antees conservation of the chemical species and energy and 
allows time dependence. The formalism also accounts for 
the interactions between the global contraction or expansion 



(with velocity v) of the star and the convective (or drifting) 
velocity u. In fact, the drifting velocities for the upward and 
downward motion (v u and t>d) do not even have to be equal. 

But it is only a formalism and requires a model for all 
horizontal motions/exchanges. In section 2.5 we suggested a 
model that is easy to implement and which can be used to 
compare our new formalism with previously derived theories 
of convection. However, this has a few limitations. 

First, we used a very simple model (where Am was set 
equal to zero) for the net mass transfer M between the up- 
ward and downward moving fluids. This eliminates the main 
non-local terms involving Am in equation (48) and hence dis- 
ables the effects of convective overshooting. To account for 
overshooting in a self-consistent way requires an a priori 
physical model for the exchange term M. However, one can 
investigate the effects of overshooting within our framework 
by using a necessarily somewhat ad hoc prescription for the 
convective velocity which allows a finite convective speed be- 
yond the formally convective region according to the Ledoux 
criterion. 

Second, we chose reversible processes for the exchange 
of momentum and the mixing of chemical species. Irre- 
versible processes would dissipate part (if not all) of the 
work done as heat. The cooling terms W conv and W{j would 
then have lower (possibly zero) values. For example, the 
Reynolds numbers in convective regions are so large that 
there is almost certainly a large amount of turbulence. To 
account properly for turbulent dissipation in the derivation 
of our model would require the inclusion of a finite viscos- 
ity introduced by the turbulent cascade. This is not easily 
done. However, one can investigate the possible outcome of 
dissipation by switching the work terms on and off. In sec- 
tion 5.3.1 (Figure 5), we investigated the effect of setting 
Wconv = and found differences of up to 20%. Furthermore, 
setting W{j = would most likely only have a small effect 
because W{j is already small. 

6.1 Stationary cores 

As a first illustration of our two-stream formalism we com- 
puted stationary convective cores. However, the convective 
core during the carbon flash is growing very rapidly. In- 
deed, the time dependence does matter at least in the en- 
ergy equation, as was shown in section 5.1. Our stationary 
cores may therefore only be very rough approximations of 
growing convective cores. Moreover, in the present study we 
used arbitrary abundances for the Urea nuclei. Most of the 
Urea nuclei are by-products of carbon burning. To compute 
their abundances in a self-consistent manner involves quite 
an extended network, as was shown by Iben (1978a). 

Despite these limitations, the results already help to 
shed some light on the convective Urea process and prove 
useful in calibrating numerical aspects in the implementa- 
tion of the method. For example, the terms AR, Ae, and 
AL, as well as the non-local term in equation (48), are found 
to be second order terms in our computations of stationary 
Urea cores. Furthermore, the chemical composition differ- 
ences between the two streams are generally small, except 
possibly at the outer edge of a convective core with a high 
Urea abundance. This may provide some justification for the 
use of the one-stream approximation to describe the convec- 
tive Urea process in a stellar evolution code. 
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6.2 Future work 

The next step is to implement the two-stream model (or 
a more simple, but appropriate approximation) in a stellar 
evolution code. Coupled with a suitable nuclear reaction net- 
work, this will allow us to follow the evolution of the core in 
a fully self-consistent manner to the runaway phase in a SN 
la and to determine the physical and chemical conditions in 
the core at the time of the explosion. This will provide the 
ignition conditions, the thermodynamic properties and the 
location of the ignition point for explosion calculations (e.g. 
Hillebrandt et al. 2004). Moreover, with this tool, we shall 
be able to systematically address the dependence of the ig- 
nition conditions on the overall metallicity, the initial C/O 
ratio, the white dwarf accretion rate, the initial mass of the 
white dwarf etc. 

In this context, we note the importance of the neutron 
excess at the time of the explosion. Timmes, Brown & Tru- 
ran (2003) have recently emphasised the role of the initial 
22 Ne abundance, which they argued was determined by the 
metallicity, introducing a metallicity dependence for SNe la. 
However, the neutron excess itself is affected by the elec- 
tron captures and emissions in the simmering phase preced- 
ing the nuclear runaway. If we take as indicative the Urea 
abundances derived by Iben (1982) (table 4) at the end of 
his computations, we can compute the number of additional 
neutrons introduced by the Urea reactions alone on Ne iso- 
topes, nu = 3F( 23 Ne) + 5Y( 25 Ne) = 3.6 x 10~ 3 . The cor- 
responding number of additional neutrons caused by 22 Ne 
resulting from solar abundances in the white dwarf progeni- 
tor is n = 2Y Q ( 22 Ne) = 2.5 x 10~ 3 . Hence the effect of the 
Urea isotopes of Ne on the neutron excess can be even larger 
than that of fossil 22 Ne. This provides another illustration 
for the importance of a proper treatment of the Urea process 
for answering some of the fundamental, unsolved questions 
concerning SNe la. 



7 SUMMARY AND CONCLUSIONS 

We have derived a two-stream formalism which carefully ad- 
dresses the energy and chemical budgets in a convective re- 
gion. In addition, it allows time dependence and describes 
the interaction of convection with the general motions of 
the star. We illustrated this formalism with a simple model 
and compared the resulting theory of convection to existing 
theories. We also derived a one-stream limit approximation 
which will be easy to implement in a stellar evolution code 
as an extension of classical MLT. 

We then applied this formalism to the convective Urea 
process and derived the entropy equation which has been 
central to previous discussions of the Urea process, and com- 
puted the convective velocity. We showed that the net heat- 
ing effect of the Urea process strongly depends on the state 
of mixing of the convective core, for which we provide an 
estimate. Urea reactions generally tend to reduce the effects 
of buoyancy. More generally we show that, in degenerate 
matter, reactions that change the number of electrons have 
a direct influence on the convective velocity. 

As an illustration of our model we computed stationary 
convective cores. These computations show that, even for a 
very small Urea fraction, convective velocities are strongly 



modified compared to the case without Urea nuclei. They 
also show that convective Urea cores are unlikely to be in a 
stationary state. Hence, time-dependent computations with 
a full nuclear reaction network are needed to provide the 
final answer to the question, "what is the influence of the 
convective Urea process on the ignition conditions in type 
la supernovae?" 



APPENDIX A: STABILITY ANALYSIS OF 
CONVECTION 

In this appendix we analyse the stationary states of con- 
vection when the evolution of the average quantities is very 
slow. We first investigate the possible available stationary 
states and then study their linear stability. 

In the following we assume that AR = Ae = AL = 
An = and neglect v and non-local terms in the difference 
equations. 



Al The stationary states 

The radiative state (u = 0) is always a stationary solution 
of our equations. The other stationary convective states are 
solutions of the system of equations (55)-(57) when u — 
has already been factored out. We now compute the solu- 
tions of this system and examine their existence and possible 
multiplicity. 



A 1.1 The cubic equation for the convective velocity 

Standard thermodynamics gives 

de d 1 ^ds dN 
or or p or or 

m ,ainT „ ainp, , dN 



Tc P ,„ „ . , <9N 
- _(V-V B ) + M '._, 



(Al) 



where s is the specific entropy, V a = ( g}°p )s,N, c p = 
(HOp.n, A* are the chemical potentials and fi' = fj, + 
T (§m)t,p- We also have 



Ah = CpAT + fi'.AN. 
Using relation (57), we get 
AT A 



^(v-v^i + z^)- 1 



(A2) 



(A3) 



where V = f^I an d H' 1 = We now identify the 

mixing length of MLT with A, so that equation (A3) becomes 
the MLT excess temperature equation. The density differ- 
ence can be written in terms of temperature and abundance 
differences if we use the thermodynamical relation 

Alnp = -(5AlnT + /x".AlnN, (A4) 

where 5 = -(|££) p>N and /x" = (|^) PjT . We obtain the 
convective velocity u by substituting this expression into 
equation (55) and using relations (57) and (A3) 

M 3 +u M 2 +u?[5(-V+V a )+//'.V N ]it+u u 2 //'.V N = 0,(A5) 
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■ X ' Ul = V P 2H P ' - in — d i n 

the dot product /x".Vn corresponds to the more familiar 
0V M where /i is mean molecular weight, 4> — ( )p,t and 

V M = . Solving this cubic equation for the real positive 
roots (if they exist) gives the convective velocity which in 
turn allows to compute all the A quantities : 

A 



o2 |-,and V N = ^ 



Note that 



AlnT = 



AlnN = 



2H P u + uq 

Vn- 

2H V 



(V - V„) and 



(A6) 
(A7) 



We now determine the number of solutions of the cubic (A5). 



A 1.2 Multiplicity of stationary states 

This cubic may have three real roots or one real root and 
two complex conjugate roots according to the sign of its 
discriminant (positive or negative). 

The number of real positive roots relies on the re- 
spective signs of <5(V — V ) — // Wn (Lcdoux criterion) 
and /x".Vn (Rayleigh- Taylor criterion). We use the rule of 
Descartes to obtain the number of solutions. 

• If /j".Vn < there is one real positive root (Rayleigh- 
Taylor instability). 

• If /x".V N > and 5(V - V a ) - m" Vn < there is no 
real positive root. 

• If £t".V N > and <5(V - V a ) - /x".V N > there is no 
positive root if the cubic discriminant is positive and there 
are two positive roots otherwise (convective case). 

In the convective case, when u ^> uo, the maximum 
convective velocity becomes 



U = UlV / 5(V-V a )-M".Vp 



(A8) 



A2.1 Radiative state (u = 0) 

In this case the matrix (A12) reads 







Sui 



X 

2H, 



V„) 

V N 



2 2Hp 
A 

•wo 




2 2H, 

o 
o 



(A13) 



Its characteristic polynomial is Pn(x) = 
— (— x) N ~ 1 P{x) where P is the cubic (A5). Hence when 
there exists a convective state (a real positive root of the 
cubic P) the radiative state must be unstable. If there is 
no real positive root to P and the discriminant is positive, 
the stability depends on whether the real part of the 
complex conjugate roots are negative. The unstable case 
corresponds to semi-convection, a state that is intrinsically 
time-dependent. In all the remaining cases (negative dis- 
criminant and no real positive root) the radiative state is 
stable. 



A 2. 2 Convective state (u > 0) 

In this case, the matrix (A12) reads 

-(u + If) 



u+uq 2H P 





-2u 
A (V- 



Va) 







2 2H P /, 

-Wi-jTM 


-u.I 



■(A14) 



When there are two available stationary states with 
different velocities u the state with the lower velocity is 
always unstable, and the state with the higher velocity is 
always stable. 

If u 3> uo the eigenvalues are real negative {— 2u, — u} and 
the convective state is stable. 



A2 Stability 
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To analyse the stability of the stationary states we need the 
time- dependent equations for the convective properties 



A dA In T X _ , , \ a i rn 

2 dt = M ^-( V - V ")~( M + M o)AlnT 



and 

A dA InN 

2 dt 



2H, 



-V N - AlnN) 



(A9) 
(A10) 

(All) 



In equation (A10) we assume that the evolution of the mean 
fluid is slow compared to the convective time scale so that 
time-derivatives of c p , [t! and p can be neglected. 

The linear stability of this system at a stationary point 
is determined by the eigenvalues of the matrix 

2 2H„ 



-2u 

£ ^(V-Va)~AlnT 
" ^V N -AlnN 



-(u + Uq) 




-ui^fl 

-u.I 



.(A12) 



where I is the identity N x N matrix (N is the total num- 
ber of species) ,Vn and A ln N are vertical A-component 
vectors, and fi" is a horizontal A-component vector. The 
system is stable when the real part of all eigenvalues is neg- 
ative. It is unstable otherwise. 
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